Excitonic effect on optical response in one-dimensional two-band Hubbard model 
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Motivated by the gigantic nonlinear optical response in the halogen-bridged Ni-compounds, the 
underlying electronic states of the compounds are examined in the one-dimensional two-band Hub- 
bard model, by studying the current-current correlation function and the charge density in the 
ground state. The dynamical density matrix renormalization group method is employed. We find 
that the low-energy peak of the correlation function consists of a single Lorentzian component for 
a parameter range appropriate to the compounds. This is due to an excitonic state induced by the 
intersite Coulomb repulsion between holes on the metal and halogen ions. This is consistent with 
the optical absorption spectra of the compounds. We suggest that the localization of holes on the 
metal ions in the ground state brings about the formation of the excitonic state. 
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One-dimensional Mott insulators such as cuprates and 
halogen-bridged Ni-compounds exhibit gigantic nonlin- 
ear optical responses, which have potential as opto- 
electronic devices. Q, |^ |^ |^ |^ It has been shown that 
this massive response arises because the photoexcited 
states across the Mott gap with even and odd parities 
are degenerate. This leads to the enhancement of the 
dipole matrix element of the third-order nonlinear sus- 
ceptibility. In the Mott insulators, both the degeneracy 
and the optical gap reflect strong electron correlation. 

The single-band Hubbard model with nearest-neighbor 
Coulomb repulsion at half-filling is a minimal model for 
the study of the optical response and photoexcited states 
in one-dimensional Mott insulators. The photoexcited 
states and optical response of the model have been exam- 
ined by employing the dynamical density matrix renor- 
malization group (DDMRG) method. |E II II i, 113, 113 
The model and the method are suitable for the clarifica- 
tion of the general electronic properties of the insulators. 
However, the real compounds such as the cuprates and 
the halogen-bridged Ni-compounds arc the insulators of 
charge transfer type, where both half-filled d-orbitals and 
empty p-orbitals for holes exist in one-dimensional net- 
works of the transition metal and ligand ions. Therefore, 
the two-band Hubbard model is relevant for the insula- 
tors. 

For the understanding of the optical response of the 
Ni-compounds with huge third-order nonlinear suscep- 
tibility, there are key electronic structure parameters in 
the two-band Hubbard model. |l|ll|ll|li| One of them 
is the charge transfer energy A, which is related to the 
optical gap. In the Ni-compounds, A is estimated to be 
comparable with the d-p hopping energy t by the angle- 
resolved photoemission spectroscopy, |l2l | and the phys- 
ical picture based on the perturbation expansion with 
respect to t/A may be violated. Another parameter is 
the intetsite Coulomb repulsion V between holes on the 
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metal and halogen ions, which causes an excitonic state. 
The presence of the excitonic states in the Ni-compounds 
has been confirmed through the optical absorption and 
photoconductivity measurements. Il3j | Since the opti- 
cal gap Wo observed in the absorption spectrum is smaller 
than the threshold energy of the photoconductivity, the 
excitation by wq is not sufficient for photocarrier gener- 
ation. This suggests that the excitonic state occurs. In 
fact, the optical absorption spectrum can be decomposed 
into a single Lorentzian component and a long tail struc- 
ture. 13 

Motivated by the optical properties of the Ni- 
compounds, we numerically calculate the current-current 
correlation function in the two-band Hubbard model. 
A numerical simulation technique applied here is the 
DDMRG method. 6, 7, 8, 9, 10, 11, 15, 16, 17J With 
changing A/t and V/t systematically, we find that the 
low-energy peak of the correlation function consists of a 
single Lorentzian component for a parameter range ap- 
propriate to the compounds. This is due to the excitonic 
state induced by V , and consistent with the optical ab- 
sorption spectra of the compounds. In addition to the 
correlation function, the charge density in the ground 
state also provides useful information to understand the 
formation of the excitonic state. Then, we calculate the 
charge density by using the DMRG method, and suggest 
that the localization of holes on metal ions in the ground 
state brings about the formation of the excitonic state. 

The two-band Hubbard model on one-dimensional 
chain is given by 

= ^ e^?^^ + ^ ii,i+i(c,[^Cj+i,o. + H.c.) 

i i.a 

+ ^C/jnjjni,| -I- y ^n^ni+i, (1) 



where the operator c| ^ creates a hole with spin a at 
site i, Ui — n^^i -I- n^^^ is the number operator, the level 
energy and the on-site Coulomb interaction Ut are de- 
fined by Ci = A and Ui = Up for odd i (halogen ion), and 
= and Ui = Ud for even i (metal ion), respectively. 
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FIG. 1: Xii^) a function of A/t and V/t. The on-site 
Coulomb repulsion is fixed to be Ud/t ~ 8 and Up/t = 0. 
Dashed lines denote single Lorentzian components. 



— t and the sign of i^^i+i is taken from the elec- 
tronic states of the Ni-compounds. We introduce the 
intersite Coulomb interaction V between holes on metal 
and neighboring halogen sites. For the charge transfer 
type insulators, A < Ud- The number of holes is equal to 
that of metal sites. For a moment, the on-site Coulomb 
interaction at halogen sites. Up, is neglected to make the 
roles of V clear, although Up may be finite in the real 
compounds. In this paper, Ud is fixed to be Ud/t = 8. 

We calculate the current-current correlation function 
defined by 
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where j = iJ2i,a^hi+ii4,cr^i+i,^ ~ H-c) is the current 
operator, |0) is the ground state with energy Eq, and 7 
is a small number. We also calculate the charge density 
in the ground state, N{i) = (0| rii |0). 

Numerical simulations to obtain x-yi'^) ^^id N{i) are 
separately carried out by the DMRG method under open 
boundary condition. The superblock is constructed by 
129 sites, i.e. 64 unit cells plus one halogen site. It should 
be noted that both ends of a one-dimensional chain are 
halogen sites, since this configuration keeps the inver- 
sion symmetry of the system after the sweep procedure 
is converged. For any parameter set, the result for 129 
site system was compared with that for smaller clusters. 



so that a boundary effect inducing localization of elec- 
trons at chain edges can be identified. Additional im- 
purity potentials on the two halogen ends are included 
by changing H with H — {V/A){ni + ul) to reduce the 
boundary effect. As will be mentioned below, the ef- 
fect produces a small hump structure at around cj = A 
for V/t > 2 in X7 even for the 129 site system. For 
calculating x-yi'^)^ the DDMRG method for the fixed 
broadening factor j/t = 0.2 is appHed. [1,0, To de- 
fine the density matrix for mixed states in the DDMRG 
method, we take three target states, i.e. the ground state 
|0), a photoexcited state j"^ |0), and the correction vector 
{LU + Ea-H + ij)~^j^ |0). For calculating N{i), the stan- 
dard DMRG method is applied. |0| For both Xj{^) ^^id 
N{i), the DMRG bases are truncated up to m = 300. 

In Fig. 1, X7(w) is plotted for various values ofA/t and 
V/t. A Lorentzian is also plotted in the Figure for V/t > 
2. It is defined by L^(w) = j'^Xii^o)/[i'^ — wo)^ + 7^] 
with uq denoting the excitation energy which provides 
the maximum x-y{^)- We focus on whether the full width 
of half maximum for a low-energy peak of x-yi^) is equal 
to that for the Lorentzian. As seen in Fig. 1, the width of 
the spectrum for A/t = 1 is qualitatively different from 
that for A/t — 2. (The spectrum for A/t = 3, which 
corresponds to a parameter for the cuprates, is similar to 
that for A/t = 2.) 

Let us start with x-yi^) for '^/t — 2 [see Figs. 1(e), 
(f), (g), and (h)]. For V/t = 0, the spectra are com- 
posed of two broad structures. One is a broad peak at 
around to/t = 1.5. This is due to the excitation pro- 
cess from the Zhan g-Ri ce singlet bound state to the up- 
per Hubbard band jig, which we call the process (A). 
The other is a shoulder at around uj/t = 3. This is 
due to the process from the unbound states to the up- 
per Hubbard band, which we call the process (B). The 
energy differences between the two structures can be 
scaled by the coupling constant of the singlet pairing, 
Jo = 4t2[(l/A) + {l/{Ud - A))]. With increasing V, the 
spectral intensity due to the process (A) increases, while 
the intensity due to the process (B) rapidly decreases. 
For V/t > 2, a main peak due to the process (A) is lo- 
cated at cjQ ~ A + y — t. The full width at half maximum 
for the main peak is almost the same as that for the single 
Lorentzian. Therefore, we find that the excitonic state 
occurs for V/t >2. It is noted that a small hump struc- 
ture seen around a; = 2 = A in Fig. 1(h) is due to the 
boundary effect. In Fig. 2, Xyi^) is replotted for various 

V with fixing A/t — 2. x-yi^a) is enhanced by increasing 

V for 2 < V/t < 3, while it tends to be suppressed with 
larger V. 

Next, we examine Xyi'^) for A/t = 1 [see Figs. 1(a), 
(b), (c), and (d)]. For V/t < 1, the intensity due to the 
process (A) increases with increasing V/t. The intensity 
due to the process (B) is almost kept. For V/t > 2, the 
main peak structure is observed near ~ A + V — t. 
However, it is not clear whether the excitonic state is 
formed or not, since the full width at half maximum 
for the main peak is larger than that for the single 
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FIG. 2: V dependence of x-,{^) for A/t = 2. 
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FIG. 3: A dependence of x(t^) for = 2. 



Lorentzian. Even though y is further increased, the fit- 
ting does not works weU. Furthermore, for Vjt > 3, 
X-y(wo) decreases with increasing V. The change of Xjii-o) 
is larger than that for A/t = 2 (not shown). An addi- 
tional peak structure due to the process (b) is also ob- 
served near uj/t = 4 in Fig. 1(c) and ui/t — 4.5 in Fig. 
1 (d) . The intensity of the additional peak increases with 
increasing V. 

As discussed in the previous two paragraphs, the exci- 
tonic state seen as a low-energy peak of X-yi^) depends 
on the magnitude of A/t. Let us examine how the ex- 
citonic state is destroyed with decreasing A/t. Since, 
in the DDMRG method, the spectral shape depends on 
the magnitude of the broadening factor 7, we remove 
the effect of the broadening by using the Lorentz trans- 
formation. After the removal of the effect, Xji'^) is 
expressed by xi^)- the actual numerical calculation, 
x(ti;) still has a finite broade ning factor (5w(< 7). The 
factor is set to be Suj = 7/2. [13 Figure 3 shows the A 
dependence of xi^) for V/t = 2. The peak intensity lin- 
early decreases with decreasing A/t. For A/t > 1.5, the 
full width at half maximum of the main peak is nearly 
twice the broadening factor, 2Suj. The width becomes 
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FIG. 4: Charge density from one halogen end ( i = ) to the 
central halogen site ( i = 65 ) for V/t — 2. (a) A/t = 3, (b) 
A/t = 2, and (c) A/t = 1. 



narrower if we choose 7 smaller than 0.2t. This indicates 
the presence of the excitonic state. However, the width 
for A/t < 1.25 is beyond 2Suj. The width increases with 
decreasing A/t for A/t < 1.25. Even for V/t > 3, the 
width for A/t ^ 1.25 is still beyond 2Slj (not shown). 

In the following, we show that the broad spectrum 
for A/t^ 1.25 is caused by delocalization of holes in the 
ground state. Figure 4 shows the charge density in the 
ground state for A/t = 3, 2 and 1 with V/t = 2. po| 
The differences of the charge density between metal and 
halogen ions for A/t = 3, 2 and 1 are 0.83, 0.71 and 0.38, 
respectively. In contrast with the difference for A/t ^ 3 
and 2, this is small for A/t = 1. The delocalization of 
holes for A/t ~ 1 is related to the destruction of the ex- 
citonic state: As seen in Fig. 3, the low-energy peak of 
x(w) for V/t = 2 consists of a single Lorentzian compo- 
nent for A/t > 1.5, while does not for A/t < 1.25. There- 
fore, the larger the charge localization is, the more stable 
the excitonic state is. In other words, holes should be 
localized on metal ions in the ground state, when the 
excitonic state is formed. 

Finally, x(^) is compared with the optical absorption 
spectrum of the Ni-compounds. For A/t > 1.5 and 
V/t = 2 in Fig. 3, x(w) is consistent with the optical 
absorption spectrum observed experimentally: x{^) f^^s 
a single Lorentzian peak followed by a long tail struc- 
ture. [3 For these parameter values, the excitation energy 
ujo which corresponds to the optical gap is estimated to be 
Wo > 2.5t - 2.5(eV), if we take t ~ l(eV). [H Compared 
with the observed optical gap [1.3(eV) ~ 2.0(eV)], loq 
may be overestimated. However, we have confirmed nu- 
merically that the finite Up extends the region where the 
excitonic state occurs. For the realistic Up and A/t^ 1, 
ujo becomes consistent with the observed optical gap. 

In summary, we have investigated the current-current 
correlation function and the charge density of one- 
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dimensional two-band Hubbard model by using the den- 
sity matrix renormalization group approach. In the two- 
band Hubbard model, the intersite Coulomb interaction 
V between holes on metal and halogen ions is a primary 
factor to produce the exciton effect on the correlation 
function. Wc found that the low-energy peak of the cor- 
relation function consists of a single Lorentzian compo- 
nent for A/t ^ 1.5 and V/t > 2 in the case that Ud/t — 8 
and Up — 0. This means the formation of the excitonic 
state. This is consistent with the optical absorption spec- 
trum observed in the Ni-compounds. By examining the 



correlation function and the charge density in the ground 
state, we suggested that the localization of holes on metal 
ions in the ground state brings about the formation of the 
excitonic state. 
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